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ABSTRACT 

Progression to malignancy requires that cells 
overcome senescence and switch to an immortal 
phenotype. Thus, exploring the genetic and epigen- 
etic changes that occur during senescence/immor- 
talization may help elucidate crucial events that lead 
to cell transformation. In the present study, we have 
globally profiled DNA methylation in relation to gene 
expression in primary, senescent and immortalized 
mouse embryonic fibroblasts. Using a high- 
resolution genome-wide mapping technique, 
followed by extensive locus-specific validation 
assays, we have identified 24 CpG islands that 
display significantly higher levels of CpG methyla- 
tion in immortalized cell lines as compared to 
primary murine fibroblasts. Several of these 
hypermethylated CpG islands are associated with 
genes involved in the MEK-ERK pathway, one of 
the most frequently disrupted pathways in cancer. 
Approximately half of the hypermethylated targets 
are developmental regulators, and bind to the re- 
pressive Polycomb group (PcG) proteins, often in 
the context of bivalent chromatin in mouse embry- 
onic stem cells. Because PcG-associated aberrant 
DNA methylation is a hallmark of several human 
malignancies, our methylation data suggest that 
epigenetic reprogramming of pluripotency genes 
may initiate cell immortalization. Consistent with 
methylome alterations, global gene expression 
analysis reveals that the vast majority of genes 
dysregulated during cell immortalization belongs 
to gene families that converge into the MEK-ERK 
pathway. Additionally, several dysregulated 



members of the MAP kinase network show con- 
comitant hypermethylation of CpG islands. 
Unlocking alternative epigenetic routes for cell 
immortalization will be paramount for understand- 
ing crucial events leading to cell transformation. 
Unlike genetic alterations, epigenetic changes are 
reversible events, and as such, can be amenable 
to pharmacological interventions, which makes 
them appealing targets for cancer therapy when 
genetic approaches prove inadequate. 

INTRODUCTION 

Cellular senescence is characterized by irreversible loss of 
cell proliferation, and is considered a powerful tumor 
suppressor mechanism to eliminate damaged cells prone 
to transformation (1,2). Cellular senescence has been 
shown to function as a barrier to somatic cell 
reprogramming (3), and to play a role in organismal 
aging (4). Features of senescence are often detected in 
pre-malignant lesions of mice and humans, but not in 
their respective invasive/malignant tumors. This suggests 
a role for senescent cells in the initiation and progression 
of tumorigenesis (2,5,6). For example, senescent cells are 
known to give rise to neighboring cancer cells through 
paracrine signaling with their microenvironment (7-9). 
Alternatively, cells may bypass or escape senescence and 
re-enter the cell cycle to become transformed following the 
inactivation of key components of the senescence- 
promoting pathways (2,4,10). Thus, investigating the 
underlying mechanisms of senescence bypass and immor- 
talization may help elucidate tumor initiation and 
progression. 

Despite obvious differences between human and mouse, 
key aspects of cellular senescence in the two species are 
well-conserved during evolution (11,12). Mouse models 
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have been shown to provide mechanistic insights into the 
chain of events that leads to senescence bypass and im- 
mortalization of relevance to human carcinogenesis 
(2,13-15). Mouse embryonic fibroblasts (MEFs) are a 
validated model system to study senescence bypass and 
spontaneous immortalization (16,17). Primary MEFs 
have a limited growth capacity and undergo stress-induced 
senescence, presumably due to oxidative damage caused 
by standard cell culture conditions (10,18,19). Upon pro- 
longed culturing and passaging, however, MEFs can 
override the senescence block and become spontaneously 
immortalized, partly due to mutation/deletion of the 
p53/pl9 Arf pathway (16,18,20). 

The phenotypically normal Hupki (human £53 knock- 
in) mouse strain, harboring the humanized p53 gene (15), 
has extensively been used as a model system to study the 
underlying mechanisms of senescence bypass and immor- 
talization (13,15,21,22). Using a library of several hundred 
immortalized cell lines, Whibley et al. (17) have 
demonstrated that MEFs derived from Hupki mice 
undergo senescence and become immortalized in the 
same way as their wild-type counterparts, i.e., embryonic 
fibroblasts from standard laboratory mice with a murine 
Trp53 gene. In both cases, immortalized cells typically 
harbor either a missense point mutation in the p53 gene 
leading to loss of p53 transactivation, or a biallelic 
deletion at the pl9 Ar ' locus. Surprisingly, however, more 
than half of the immortalized cell lines examined had 
neither a p53 mutation nor loss of the pl9 Arf activity 
(17). In some of these cell lines, genetic defects in tumor 
suppressor pathways, other than those affecting the 
p53/pl9 Arf axis, may contribute to escape or release from 
the senescence block into an immortal phenotype (17,23). 
Furthermore, active acquisition of epigenetic modifica- 
tions, such as DNA methylation at regulatory CpG sites 
may interfere with the expression of key senescence genes, 
and pave the way to immortalization and cell transform- 
ation (23,24). To explore novel epigene tic/genetic routes 
that may modulate senescence bypass and immortaliza- 
tion, we have globally profiled DNA methylation in 
relation to gene expression in primary, senescent and 
immortalized fibroblasts derived from both the Hupki 
and wild-type mice. Using microarray-based approaches 
together with locus/gene-specific validation assays, we 
have constructed the whole methylome and transcriptome 
of MEFs at three distinct proliferation stages, including 
primary, senescent and immortalized cell lines. 

MATERIALS AND METHODS 

Isolation and immortalization of primary MEFs 

Primary MEFs were obtained from embryos derived from 
either Hupki (Sets 1, 2 and 3) or wild-type (Sets 5 and 7) 
mice according to standard procedures (15). Briefly, 
mouse embryos were harvested in utero at 13.5 days of 
gestation. Following the removal of the head and 
internal organs, embryo carcasses were roughly minced 
and digested for five times, 5min each, with 0.25% 
trypsin at 37°C. Remaining tissue clumps were 
disaggregated by passage through a syringe equipped 



with an 18-gauge needle, and vigorous pipetting. 
Homogenous cell suspensions were then transferred to 
tissue culture flasks, and cultured in standard Dulbecco 
Modified Eagle's Medium (DMEM, Irvine Scientific, 
Santa Ana, CA), supplemented with 10% fetal calf 
serum (FCS), non-essential amino acids and 0.1 mM 
(3-mercaptoethanol. Following 1-3 days of incubation at 
37° C, confluent cultures were trypsinized and frozen, at 
2-3 x 10 6 cells per vial. To establish cell lines, 
cryopreserved embryonic fibroblasts (passage 0) were 
thawed and seeded onto T75 flasks containing DMEM 
supplemented with 10% FCS. At 80% confluence, the 
cells were trypsinized and re-seeded onto 6-well culture 
dishes or small flasks (12.5 cm 2 ) at a density of ~1 x 10 5 
cells per well/flask in DMEM plus 10% FCS. Cultures 
were maintained under standard conditions (37°C, 5% 
C0 2 in a humidified atmosphere), and passaged serially 
when reaching confluence [pass to: 1:3 or 1:4, as described 
previously (17)]. During crisis, which typically occurs by 
passage 5-6 and lasts up to 3 weeks, culture medium was 
changed twice a week. Following the crisis, over 95% of 
all cultures became spontaneously immortalized. For all 
experiments, cells were harvested at different passages 
before-, during- and after crisis, as outlined in 
Figure 1A, and genomic DNA/RNA was isolated accord- 
ing to standard procedures (17). 

Genome-wide DNA methylation profiling by 
MIRA-microarray analysis 

We used the methylated CpG (mCG) island recovery 
assay (MIRA) in combination with microarray analysis 
(25,26) to globally profile DNA methylation in primary, 
senescent and immortalized MEFs. As a pull-down assay 
for enrichment of mCG islands, the MIRA is based on the 
ability of the methyl-CpG binding 2b (MBD2b) protein to 
bind methylated-CpG dinucleotides, while this reaction is 
enhanced in the presence of MBD3L1 protein (27). The 
MIRA-enriched- and input (non-enriched) DNA fractions 
were amplified by polymerase chain reaction (PCR), and 
subsequently labeled and hybridized to the Roche 
NimbleGen mouse methylation microarrays (Roche 
NimbleGen Inc., Indianapolis, IN) to interrogate the 
entire CpG-island battery of the mouse genome (see 
Supplementary Data) (25,26). The accession number for 
DNA methylation microarray data is GSE39034. 

Single gene I locus DNA methylation analysis by 
COBRA and bisulfite sequencing 

We verified the methylation status of individual target 
loci/genes identified by MIRA-microarray analysis using 
the conventional combined bisulfite restriction analysis 
(COBRA) (28), and sodium bisulfite sequencing (29). 
The primer sets used for PCR amplification of all the 
analysed targets are listed in Supplementary Table SI 
(see Supplementary Data). 

DNA methylation analysis of repetitive DNA elements 
We have specifically determined the methylation status of 
major repetitive DNA elements, including long 
interspersed nuclear elements (LINE LI), intracisternal 
A particle long terminal repeat retrotransposons 
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Figure 1. Outline of the passaging protocol and identification of differentially methylated targets in immortalized MEFs. (A) MEFs were cultured 
in vitro as described in the text and harvested at low-passage (before crisis/primary), during crisis (crisis/senescent), and following immortalization 
(after crisis/immortalized). DNA/RNA was isolated according to standard procedures (17). (B) Aberrantly mCG islands were identified by 
MIRA-microarray analysis. Three independently established immortalized murine cell lines are shown (Set 1, Set 2 and Set 3). CpG islands were 
considered as hypermethylated when the immortalized cell line showed methylation peaks that were absent in all controls (primary MEFs). The 
minimum difference between immortalized and primary MEFs was >2-fold. For hypomethylated targets, reverse criteria were used. The numbers of 
methylation peaks, per each dataset, are shown. 



(IAP-LTR), and short interspersed nuclear elements 
(SINE Bl) (30-32), in the genome of primary, senescent 
or immortalized MEFs using a bisulfite-based sequencing 
analysis (see Supplementary Data) (33). 

Genome-wide gene expression profiling by 
Affymetrix-microarray analysis 

We used the GeneChip® Mouse Genome 430 2.0 Array 
(Affymetrix Inc., Santa Clara, CA) to establish the global 
gene expression profile of MEFs at three different stages 
(before crisis, crisis, after crisis). Currently, this micro- 
array platform is the most comprehensive mouse 
genome expression array, which enables interrogation of 
over 39 000 transcripts and variants from over 34 000 
well-characterized mouse genes (see Supplementary 
Data). The accession number for gene expression micro- 
array data is GSE39034. 



Quantitative real-time PCR 

Standard quantitative real-time-PCR (qRT-PCR) was 
used to confirm the transcriptional status of individual 
target genes identified by microarray analysis (34). The 
primer sets used for qRT-PCR are listed in 
Supplementary Table SI (see Supplementary Data). 



RESULTS 

Genome-wide detection of aberrant DNA methylation 
following senescence bypass and immortalization 

We have used the MIRA in combination with microarray 
analysis (25,26) to detect changes in the global methyla- 
tion pattern of murine cells undergoing immortalization. 
MEFs were cultured in vitro as described in 'Materials and 
Methods' section (outlined in Figure 1 A), and harvested at 
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different passages (i.e., before crisis (be): primary, during 
crisis (c): senescent and after crisis (ac): immortalized). 
Following pull-down of mCG islands by MIRA, the 
enriched- and input (non-enriched) DNA fractions were 
labeled, mixed and hybridized to the mouse CpG island 
plus promoter tiling arrays (Roche NimbleGen Inc.). 
Using standard algorithms for peak calling (see 
Supplementary Data), we have identified 24 CpG islands 
that were aberrantly methylated in three independent sets 
of immortalized MEF cell lines, but were methylation-free 
in the correspondent low-passage primary cells (Table 1). 
The identified CpG islands mapped to intergenic regions 
(21%), within gene bodies (38%), at promoter regions 
(29%) and near the 3'-end of genes (12%). Using a com- 
bination of the Ingenuity Pathway Analysis® (IP A®: v 9.0) 
and the Database for Annotation, Visualization and 
Integrated Discovery (DAVID: v 6.7) (35), we obtained 
gene ontology (GO) information for the annotated genes 
identified as hypermethylated targets in immortalized 
MEFs. Functional annotation analysis revealed that 
gene targets involved in intracellular signaling cascade 
{Rem 2, Rasdl, Wnk4, Adrald, Mapkll, Plcd3; 
P< 0.007), small GTPase-mediated signal transduction 
(Rem2, Rasdl, Mapkll; P < 0.04), as well as neuropeptide 
signaling pathway (Eltdl, Lphn3; P<0.09) were particu- 
larly enriched. Of significance, several genes involved in 
the MEK-ERK pathway, one of the most frequently dis- 
rupted pathways in cancer (36-39), showed progressive 
DNA hypermethylation as the cells bypassed senescence 
and became immortalized. Using reverse criteria for data 
analysis, we detected no significant loss of CpG island 



methylation across the three independently established 
immortalized cell lines relative to counterpart primary 
cells. 

Verification of aberrantly methylated targets by 
COBRA and bisulfite sequencing 

To validate the DNA methylation microarray results, we 
have randomly selected and analysed several of the targets 
identified by MIRA-microarray analysis using both the 
COBRA and bisulfite DNA sequencing analysis (28,29). 
Representative results of the COBRA assay are shown in 
Figures 2C and 3C and Supplementary Figures S1-S4. 
Here, we have inspected the methylation status of CpG 
islands associated with the Plcd3 gene (target #2), Phox2a 
(target #6), Mapkll (target #9), Cldn5 (target #11) and 
Rasdl (target #18) genes. As shown in Table 1, these 
CpG-rich loci displayed enrichment ratios, by MIRA- 
microarray analysis, ranging from 5.99 to 3.48, respect- 
ively (on a linear scale, immortalized versus low-passage 
fibroblasts). COBRA analysis of the above-specified 
targets confirms the microarray results, indicating that 
novel epigenetic marks, i.e., CpG methylation, are estab- 
lished during cell immortalization (Figures 2A and 3A and 
Supplementary Figures S1A-S3A). Of note, these methy- 
lation changes were consistent in both the Hupki and 
wild-type immortalized MEFs (Supplementary Figures 
S3 and S4). 

To obtain a detailed methylation map for each CpG 
within the CpG islands, we also performed bisulfite 
sequencing on the targets analysed by COBRA. The 
sequencing results of the Phox2a-, Mapkll-, Rasdl- and 



Table 1. Hypermethylated targets detected in immortalized versus primary fibroblasts 



Target no. 


Gene symbol 


Description 


CGI location 


Position relative 
to gene 


Strand 


Enrichment 
ratio 


1 






chrl9:59543322-59543961 






6.40 


2 


Plcd3 


phospholipase C, delta 3 


chr 1 1 : 1 0294 1383-1 02942047 


Intragenic 




5.99 


3 


Wnk4 


WNK lysine deficient protein kinase 4 


chrll:101126459-101126913 


Intragenic 


+ 


5.09 


4 


Adrald 


adrenergic receptor, alpha Id 


chr2:131371751-131372295 


Downstream 




5.01 


5 






chrl3:56854334-56854768 






4.85 


6 


Phox2a 


paired-like homeobox 2a 


chr7: 108969059-108969729 


Intragenic 


+ 


4.80 


7 


Rem2 


Rad- and gem-related GTP-binding 


chrl4:55098125-55098574 


Downstream 


+ 


4.78 






protein 2 










8 






chr7:77514032-77515167 






4.71 


9 


Mapkll 


mitogen-activated protein kinase 11 


chrl5:88976796-88977465 


Intragenic 




4.67 


10 


Sh2d5 


SH2 domain containing 5 


chr4: 137806496- 137806960 


Upstream 


+ 


4.25 


11 


Cldn5 


claudin 5 


chrl6:18777229-18778053 


Intragenic 


+ 


4.20 


12 


Hoxc9 


homeobox C9 


chrl5:102807421-102808173 


Upstream 


+ 


4.20 


13 


Ntn3 


netrin 3 


chrl7:24345612-24346041 


Upstream 




4.19 


14 


Fbrsll 




chr5:l 10815430-1 10815859 


Intragenic 




4.16 


15 


Ddn 


dendrin 


chrl5:98637075-98637514 


Intragenic 




3.85 


16 


Amh 


anti-Mullerian hormone 


chrl0:80269858-80270292 


Downstream 


+ 


3.67 


17 


Eltdl 


EGF, latrophilin seven transmembrane 


chr 3 : 1 5 1 1 00756- 1 5 1 1 0 1 200 


Upstream 


+ 


3.57 






domain containing 1 










18 


Rasdl 


RAS, dexamethasone-induced 1 


chrl 1:59777219-59777763 


Intragenic 3 




3.48 


19 


Adrald 


adrenergic receptor, alpha Id 


chr2:131386864-131387318 


Upstream 




3.46 


20 


Lphn3 


latrophilin 3 


chr5:81758777-81759233 


Intragenic 


+ 


3.38 


21 


Akap5 


A kinase (PRKA) anchor protein 5 


chrl2:77425710-77426441 


Upstream 


+ 


3.19 


22 


CcdclU 


coiled-coil domain containing 114 


chr7:53179579-53180028 


Upstream 


+ 


3.06 


23 






chrl 1:21903730-21904098 






2.96 


24 






chrl5:102795737-102796201 






2.86 



a Downstream of Med9 gene. 
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Figure 2. Hypermethylated CpG island associated with the Phox2a locus. (A) Visualization of the NimbleGen array data for immortalized (in red) 
and primary MEFs (in green). Three independently established immortalized murine cell lines are shown (Set 1, Set 2 and Set 3). At this level of 
resolution, each peak corresponds to a CpG island. The methylation signal, obtained by the NimbleScan software, is plotted along the chromosome 
as a P-value on the y-axis (starting from 0 when P-value is 1). The P-value is derived from the Kolmogorov-Smirnov test comparing the log2 
enrichment ratios between MIR A and input within a 750 bp window. (B) Snapshot of the UCSC Genome Browser indicating the location of the 
Phox2a gene. CpG islands are shown as green boxes. The locations of the H3K4me3, CTCF and Pol2-binding sites, are shown for low-passage 
MEFs, based on ENCODE/LICR published databases (http://genome.ucsc.edu/ENCODE/downloadsMouse.html). H3K4me3 and Pol2 signals are 
strong indicators of active promoters, whereas CTCF-binding sites are considered as a mark for potential insulator elements. (C) Genomic DNA 
from primary (lbc, 2bc, 3bc), senescent (2c, 3c) and immortalized (lac, 2ac, 3ac) MEFs was treated with sodium bisulfite, and the CpG island within 
the Phox2a gene was amplified with gene-specific primers, and subjected to COBRA. Digested fragments on the gel are indicative of methylated 
BstUI restriction sites (5 -CGCG) within the CpG island. Mouse genomic DNA was methylated in vitro with the SssI methyltransferase, and served 
as positive control (Pos). (+) and (— ) represent the presence and absence, respectively, of the BstUI restriction enzyme in the reaction mix. 
M = 100 bp ladder DNA marker. Given the limited amount of DNA available from sample lc, bisulfite sequencing was performed only. (D) The 
extent of CpG methylation within the Phox2a CpG island was determined by sodium bisulfite sequencing in primary (lbc, 2bc, 3bc; left panel), 
senescent (lc, 2c, 3c; middle panel) and immortalized (lac, 2ac, 3ac; right panel) MEFs. The sequencing results of up to 10 independent clones, and 
the respective percentage of methylation per sample are shown. Data were analysed using the web-based software, Qu antification tool for 
Microarray Analysis (QUMA) (http://quma.cdb.riken.jp/). Open and closed circles represent unmethylated and mCG dinucleotides, respectively. 



CMz5-associated CpG islands are shown in Figures 2D 
and 3D, and Supplementary Figures SID and S2D. In 
all cases, we confirmed DNA methylation differences (up 
to 6-fold) between immortalized and primary MEFs, mir- 
roring the differences observed in the microarray and 
COBRA analyses (Table 1). As shown in Figures 2D 
and 3D and Supplementary Figures SID and S2D, 
aberrant methylation at the CpG islands is established 
mostly during the crisis stage, suggesting that DNA 
methylation of crucial genes involved in the maintenance 
of senescence may drive cells to the immortal phenotype. 



We did not detect any gene-specific CpG island that was 
hypomethylated across all the three immortalized cell 
lines. However, when less stringent criteria were used in 
the analysis (a methylation peak present in all primary 
cells but absent in two out of the three immortalized cell 
lines), we found several hypomethylated CpG islands 
common to two immortalized cell lines. Representative 
examples of hypomethylated targets are shown in 
Figure 4, which combines bisulfite sequencing results for 
the Neatl gene (Figure 4A) and COBRA data for an 
unknown RIKEN cDNA 4833420G17 gene (Figure 4B) 
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Figure 3. Hypermethylated CpG island associated with the Mapkll locus. (A) Visualization of the NimbleGen array data for immortalized (in red) 
and primary (in green) MEFs. Three independently established immortalized murine cell lines are shown (Set 1, Set 2 and Set 3). The methylation 
signal is plotted as a P-value on the y-axis, and was calculated as described in legend for Figure 2. (B) Snapshot of the UCSC Genome Browser 
indicating the location of the Mapkll gene. The corresponding CpG island is shown as a green box together with the overlapping binding site for 
CTCF, and potential transcriptional initiation sites. (C) COBRA analysis of the mCG island associated with the Mapkll gene in primary (lbc, 2bc, 
3bc), senescent (2c, 3c) and immortalized (lac, 2ac, 3ac) MEFs. Digested fragments on the gel are indicative of methylated BstUI restriction sites 
^-CGCG) within the CpG island (see, legend for Figure 2). (D) The extent of CpG methylation within the Mapkll CpG island was determined by 
sodium bisulfite sequencing in primary (lbc, 2bc, 3bc; left), senescent (lc, 2c, 3c; middle) and immortalized (lac, 2ac, 3ac; right) MEFs. 
The sequencing results of several independent clones, and the respective percentage of methylation per sample are shown. 



and the Hoxa3 gene (Figure 4C, bottom). As illustrated in 
Figure 4C, loss of methylation may occur at entire gene 
clusters, e.g., the Hoxa homeobox gene cluster that is 
known to be aberrantly methylated in a variety of 
human cancers (25,40,41). 

Using a sodium bisulfite-based sequencing analysis (33), 
we have also investigated the methylation status of major 
repetitive DNA elements, including LINE LI, IAP-LTR 
and SINE Bl, which constitute 18.78, 3.13 and 2.66%, 
respectively, of the mouse genome (30-32). As shown in 
Figure 5, DNA methylation at these repeat elements, in 
particular IAP-LTR and SINE Bl elements, was reduced 
in immortalized cell lines relative to primary MEFs, which 
indicates that genome-wide DNA methylation decreases 
with in vitro passaging. 

Gene expression profiling of primary, senescent 
and immortalized MEFs 

To identify novel genes that are differentially expressed 
during immortalization, we have profiled the whole 



transcriptome of primary, senescent and immortalized 
MEFs using an Affymetrix-microarray platform 
(Affymetrix Inc.). Following standard bioinformatics 
analysis, we identified a large number of dysregulated 
genes across three independently established immortalized 
cell lines (Figure 6A). When we compared the expression 
profile of senescent versus primary fibroblasts (based on a 
log2 ratio of one = 2-fold difference), we detected 95 
differentially expressed targets/probes, of which 34 were 
overexpressed (up to a log2 ratio of +5.54), and 61 were 
down-regulated (with a log2 ratio as low as —3.51) in all 
the three datasets. A much higher number of differentially 
expressed targets/probes (n = 313) was detected when 
immortalized cells were compared to low-passage 
primary fibroblasts (Figure 6B and Supplementary Table 
S2). Of the 313 identified targets, 80 were up-regulated 
(up to a log2 ratio of +6.07) and 233 were down-regulated 
(with a log2 ratio as low as —3.66). Interestingly, 50 of 
these probes/targets were common in the senescent and 
immortalized phenotype, suggesting a crucial role for 
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Figure 4. Hypomethylated CpG islands associated with the Neatl, RIKEN cDNA 4833420G17 and Hoxa3 locus. (A) Loss of CpG methylation 
within the Neatl CpG island was determined by sodium bisulfite sequencing in primary (lbc, 2bc, 3bc; left panel), senescent (lc, 2c, 3c; middle panel) 
and immortalized (lac, 2ac, 3ac; right panel) MEFs. The sequencing results of several independent clones and the respective percentage of methy- 
lation per sample are shown. (B) COBRA analysis of the hypomethylated CpG island associated with an unknown RIKEN cDNA 4833420G17 gene 
in primary (lbc, 2bc, 3bc), senescent (2c, 3c) and immortalized (lac, 2ac, 3ac) MEFs. Digested fragments on the gel are indicative of methylated TaqI 
restriction sites (5'-TCGA) within the CpG island (see, legend for Figure 2). (C) Upper panel: hypomethylated CpG islands associated with the Hoxa 
cluster. Visualization of the NimbleGen array data for immortalized (ac, in red) and primary (be, in green) MEFs. The methylation signal is plotted 
as a P-value on the y-axis. The location of the hypomethylated CpG island analysed by COBRA is shown as a snapshot of the UCSC Genome 
Browser. Lower panel: COBRA analysis of the hypomethylated CpG island associated with the Hoxa3 locus in primary (lbc, 2bc, 3bc), senescent 
(2c, 3c) and immortalized (lac, 2ac, 3ac) MEFs. Digested fragments on the gel are indicative of methylated Bst UI restriction sites (5'-CGCG) within 
the CpG island (see, legend for Figure 2). 



these genes in circumventing senescence and initiating un- 
controlled cell proliferation (Figure 6B). 

We have also investigated the functional relationships 
between genes that were differentially regulated in 
immortalized versus primary fibroblasts. Based on the 
IPA®, we found that these genes were mostly involved in 
cell signaling networks, cell cycle and cell morphology, 
embryonic and tissue development (Figure 6D) and were 
overrepresented compared to the senescent phenotype 
(Figure 6C). In agreement with the methylation data, we 
found that several pathways implicated in human cancer, 
such as the MEK-ERK, PI3K/Akt, Nfkb and Vegf 
pathway (36-39) were mostly affected upon immortaliza- 
tion (Figures 7 A and 8). These findings are in agreement 
with those of a recent report in which major pathways 
involved in controlling senescence/immortalization were 
identified (42). 

Correlation between DNA methylation and gene 
expression during senescence bypass and immortalization 

To investigate the impact of DNA methylation on gene 
expression during senescence bypass and cell 



immortalization, we have individually quantified the tran- 
scription level of several of the hypermethylated genes 
identified by our DNA methylation analysis using 
standard qRT-PCR. The analysed targets are known to 
interact with key components of the MAP kinase 
pathways (Figure 7A), and are potentially relevant for 
the establishment of a malignant phenotype. Figure 7B 
shows the mean normalized expression levels of the 
Cldn5, Hoxc9, Wnk4, Mapkll, Rasdl, Adrald and 
Akap5 genes in three independently immortalized cell 
lines (lac, 2ac, 3ac) relative to primary cells (lbc, 2bc, 
3bc). With the exception of the Akap5 gene, which is 
down-regulated, all the other targets show steady 
increase in transcription levels upon immortalization. 

We have also selected several of the up-regulated and 
down-regulated genes detected by Affymetrix-chip 
analysis, and validated their relative expression by 
standard qRT-PCR. Representative qRT-PCR results 
are shown in Figure 7C. In line with microarray data 
(Supplementary Figure S5 and Supplementary Table S2), 
we confirmed that the Plfl and Foxgl genes were 
up-regulated, whereas the HoxdlO, Hoxdl3 and Igf2 
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Figure 5. Methylation profiling in repetitive DNA elements in the mouse genome during immortalization. (A) Bisulfite sequencing of LINE LI, IAP-LTR 
and SINE Bl elements. Two independently established immortalized murine cell lines are shown (Set 5 and Set 7). be = before crisis; ac = after crisis; 
5 and 7 designate Set 5 and Set 7, respectively. Values in parentheses indicate percentages of mCG. (B) Quantification of mCG in major repetitive DNA 
elements in immortalized versus primary MEFs. Average results from two independently established immortalized murine cell lines are shown (Set 5 and 
Set 7). Fisher's exact test was used to calculate the statistical significance of difference in mCpGs between immortalized and primary MEFs. 



genes were down-regulated in immortalized MEFs 
(Figure 7C). Of significance, some of these genes encode 
proteins, which are relevant for MAPK cascades 
(Figures 7 A and 8). 

DISCUSSION 

Bypass from senescence is widely recognized as the first 
step toward immortalization, which, in turn, may drive 
cells to malignancy (43). To gain insights into the epigen- 
etic and genetic mechanisms involved in the senescence/ 
immortalization switch, we have analysed the whole 
methylome and transcriptome of primary, senescent and 
immortalized MEFs, a well-established in vitro model of 
cell immortalization (16,17). Using a high-resolution 
genome-wide mapping technique (MIRA-chip) (25,26), 
followed by extensive single gene/locus-specific validation 
assays (COBRA and bisulfite sequencing) (28,29), we have 



identified 24 CpG islands that display aberrant DNA 
methylation in three independently established 
immortalized MEF cell lines as compared to counterpart 
low-passage primary cells (Table 1). Furthermore, we 
observed significant loss of CpG methylation at major 
repetitive DNA elements (Figure 5), a genome- wide 
epigenetic event often associated with genetic instability 
and cancer (44-46). 

In immortalized murine fibroblasts, de novo DNA 
methylation marks preferentially genes that are involved 
in the development and function of the nervous system 
(Phox2a, target #6; Rem2, target #7; Cldn5, target #11; 
Ntn3, target #13; Amh, target #16; Adrald, target 
#19; Akap5, target #21; IPA®, P = 9.75E-04 - 4.30E-02; 
Table 1). These and other developmental regulators (47% 
of the identified genes) are targets of the Polycomb repres- 
sive complexes, and often associate with bivalent modifi- 
cations (H3K4me3 active marks and H3K27me3 inactive 
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marks) in murine embryonic stem (ES) cells, as deduced 
by comparison of our data to published databases (47-49). 
Polycomb group (PcG) proteins are critical for the main- 
tenance of stem cell identity and cell differentiation 
(47,50,51), and PcG-associated genes in ES cells are 
often the target of aberrant CpG methylation in cancer, 
thus, supporting the 'stem cell origin' model of transform- 
ation (52-54). The DNA methylation/PcG connection has 
been consistently observed in the epigenome of several 
human malignancies (55,56), and our findings suggest 
that epigenetic reprogramming of pluripotency genes 
may take place at early pre-neoplastic stages, once the 
cells acquire an infinite life span. Interestingly, GO 
analysis of the annotated hypermethylated targets 
(n = 19), identified in immortalized MEFs, confirmed a 
trend of enrichment for genes involved in intracellular 



signaling cascades, such as the MEK-ERK pathway, 
one of the most frequently altered genetic pathways in 
human cancer (Figures 7A and 8) (36-39). This finding 
reinforces the notion that alternative epigenetic mechan- 
isms may contribute to the disruption of genetic pathways 
that are crucial for cell immortalization and growth, and 
tumor development (2,3,17,18,42). 

To examine the potential impact of DNA methylation 
on gene expression during senescence bypass and immor- 
talization, we have also constructed the whole transcrip- 
tome of primary, senescent and immortalized murine 
fibroblasts. A large number of differentially expressed 
targets/probes (n = 313) was identified in independently 
established immortalized MEF cell lines when compared 
to low-passage primary fibroblasts (Figure 6B and 
Supplementary Table S2). Apart from few known 
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Figure 7. Gene network analysis and relative quantification of gene expression by standard qRT-PCR. (A) Upon immortalization, genetic and epi- 
genetic events may contribute to dysregulation of MAP kinase pathways. Upstream and downstream components of the p38-MAPK and MEK-ERK 
cascade are shown in pink. Network interactions were deduced from IP A® analysis and are illustrated in detail in Figure 8. Continuous line, direct 
interaction; broken line, indirect interaction. Red and green molecules represent up-regulated and down-regulated targets, respectively. Yellow node 
color is indicative of hypermethylation. Up-white-arrows represent up-regulation and down-white-arrows represent down-regulation. (B) Transcription 
levels of the Cldn5, Hoxc9, Wnk4, Mapkll, Rasdl, Adrald and Akap5 genes following immortalization. Bars represent the mean normalized expression 
(±SD) of three independent immortalized cells lines (ac: lac & 2ac & 3ac; each in triplicate) relative to primary cells (be: lbc & 2bc & 3bc; each in 
triplicate). Data analysis was performed using the 2~ AACt method (34) with Gapdh as reference gene. (C) The expression status of individual target 
genes identified by Affymetrix-microarray analysis was confirmed by standard qRT-PCR. Two up-regulated (Plfl and Foxgl) and three down-regulated 
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as calibrators within the respective group (fold ratio =1). Data are expressed as mean ± SD of three replicates. 
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markers of senescence [i.e., Igfbp6, Ahr, Collal, Aldhla2; 
reviewed in (42)], several of our identified transcripts are 
novel, and represent starting points for future research on 
transcriptional pathways that may regulate the acquisition 
of an infinite life span in somatic cells. We have further 
classified the function of the differentially expressed genes 
(crisis and after crisis versus before crisis phenotypes) by 
GO analysis. According to DAVID, a larger number of 
dysfunctional gene clusters was identified in the immortal 
phenotype, including genes involved in signaling, develop- 
ment, morphogenesis and neurogenesis (Figure 6C and 
D). The vast majority of these dysregulated genes 
converge into specific pathways, such as MEK/ERK, 
PI3K/Akt, NFkB and Vegf that are often impaired 
during tumorigenesis (Figure 8) (36-39). In addition, out 
of the 244 annotated DAVID IDs (from the 313 tran- 
scripts differentially expressed in immortalized MEFs), 



41% (n = 101) are normally expressed in the brain, and 
play a role in neural development. 

We have also confirmed gene expression changes in 
several of the hypermethylated targets identified by our 
DNA methylation analysis using standard qRT-PCR. 
Apart from Akap5 and Hoxc9, which are aberrantly 
methylated at their promoter regions, the majority of the 
analysed targets contain intragenic CpG islands that 
become hypermethylated upon immortalization. These 
latter genes are all up-regulated in immortalized versus 
low-passage primary cells (Figure 7B). Our findings are 
in agreement with several reports indicating that 
intragenic methylation often correlates with gene up- 
regulation [reviewed in (57-61)]. Along with Cldn5, 
Wnk4, Mapkll and Rasdl, several other targets are 
methylated within gene bodies (38%) rather than at 
canonical 5'-end promoters (29%) (Table 1), a trend 




Figure 8. Functionally related pathways of aberrantly methylated and differentially expressed genes in immortalized MEFs. Merging Networks were 
generated by Ingenuity Pathway Analysis® (IP A®: v 9.0), using the gene lists in Table 1 and Supplementary Table S2 (ac versus be). The convergence 
of genes involved in the MEK-ERK pathway, one of the most frequently disrupted pathways in human cancer (36-39), is highly notable. Red and 
green nodes represent up-regulated and down-regulated genes, respectively. Yellow nodes are indicative of hypermethylated targets and white nodes 
show molecules that are not included in the datasets but interact with other components of the networks. Canonical pathways (CP: Molecular 
Mechanisms of Cancer) are outlined in blue. 
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consistently reported for different human epigenomes 
(62,63). Although the exact functional role of these 
intragenic CpG-rich elements remains to be determined, 
the general consent is that these elements constitute short 
independent transcriptional units and/or regulate alterna- 
tive isoforms (62,63). A recent study has demonstrated a 
crucial role for the DNA-binding protein, CCCTC 
-binding factor (CTCF), in regulating alternative 
exon splicing by inhibiting Pol2 elongation dynamics 
(64). CTCF binding to the DNA occurs mostly 
intragenically (40-45%), and is inhibited by CpG methy- 
lation (64). In the present study, we have observed that 
about 70% of the methylated intragenic CpG islands are 
located at or adjacent to (within 1 kb) potential consensus 
elements for CTCF, when compared to published 
databases (Figures 2B and 3B, and Supplementary 
Figures SI and S3) (http://genome.ucsc.edu/ENCODE/ 
downloadsMouse.html). It is plausible that de novo methy- 
lation established during cell immortalization may impede 
the binding of the CTCF protein at these intragenic sites, 
thereby, resuming RNA polymerase extension and/or 
splicing. Alternatively, methylation-dependent CTCF 
release within gene bodies can predispose the DNA to 
an open chromatin configuration, which may lead to 
unusual gene activation, uncontrolled cell proliferation 
(i.e., immortalization), and cancer development (65,66). 

A large body of evidence exists indicating that the Rasj 
Raf/ERK and the p38 MAP kinase pathways are two of 
the most relevant pathways involved in the activation 
of senescence as a fail-safe mechanism against oncogenic 
transformation and other damage in primary cells 
(42,67-69). By marked contrast, deregulated activation 
of the MEK/MAPK cascade can result in forced mito- 
genesis and transformation in immortal cells (67). These 
two totally opposite outcomes, cell cycle arrest and induc- 
tion of proliferation, imposed by identical signal- 
transduction pathways are mostly cell context-dependent 
and dictated by the strength of the stimuli and the integ- 
rity of the senescence program (67). Interestingly, 
high-passage MEFs (PD 13), cultured according to the 
3T3 protocol (as in this study), displayed elevated levels 
of activated p38 MAPK (p38-P), compared to early passage 
MEF cells (PD 3), although apparently immortalized 
MEFs (PD 35) proliferate independently of the status of 
p38 MAPK (70). In our study, although we did not directly 
analyse the activity (i.e. phosphorylation) of the kinases 
involved in the MAPK signaling cascades, all the indica- 
tions are that these cascades are indeed activated during 
immortalization. Of note, we have shown that proliferin 1 
(Plfl) a well-known component of the MAPK pathway is 
greatly up-regulated in senescent and immortal fibroblasts 
(up to 67-fold on linear scale; Supplementary Table S2). 
Plfl modulates the activation of the MEK/ERK pathway 
and is transcriptionally regulated by fi-catenin, which 
accumulates following activation of p38 MAPK , via repres- 
sion of GSK3/3 (Figure 7A) (71,72). 

In conclusion, our global profiling of the methylome 
and transcriptome of primary, senescent and immortalized 
murine fibroblast cells shows abrogations of epigenetic 
and genetic pathways during cell immortalization, 
which closely resemble those found in human cancer 



(36,37,39,57-59). Despite species differences between 
human and mouse, identification of novel biomarkers of 
senescence bypass and immortalization in mouse models 
may help elucidate the underlying mechanisms of human 
carcinogenesis. If validated in vivo in humans, these bio- 
logical markers can be used to predict cancer risk in early 
stage lesions, and help devise future strategies for cancer 
therapies, e.g., re-activation of key signaling effectors of 
senescence in tumor cells. 
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